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ABSTRACT 

X-ray observations of solar flares routinely reveal an impulsive high-energy 
and a gradual low-energy emission component, whose relationship is one of the 
key issues of solar flare study. The gradual and impulsive emission components 
are believed to be associated with, respectively, the thermal and nonthermal 
components identified in spectral fitting. In this paper, a prominent ~ 50 second 
hard X-ray (HXR) pulse of a simple GOES class C7.5 flare on 20 February 2002 
is used to study the association between high energy, non-thermal and impulsive 
evolution, and low energy, thermal and gradual evolution. We use regularized 
methods to obtain time derivatives of photon fluxes to quantify the time evolution 
as a function of photon energy, obtaining a break energy between impulsive and 
gradual behavior. These break energies are consistent with a constant value of 
~ 11 keV in agreement with those found spectroscopically between thermal and 
non-thermal components, but the relative errors of the former are greater than 
15% and much greater than the a few percent errors found from the spectral 
fitting. These errors only weakly depend on assuming an underlying spectral 
model for the photons, pointing to the current data being inadequate to reduce 
the uncertainties rather than there being a problem associated with an assumed 
model. The time derivative method is used to test for the presence of a 'pivot 
energy' in this flare. Although these pivot energies are marginally consistent 
with a constant value of ~ 9 keV, its values in the HXR rise phase appear to 
be lower than those in the decay phase. Assuming that electrons producing 
the high-energy component have a power law distribution and are accelerated 
from relatively hot regions of a background plasma responsible for the observed 
thermal component, a low limit is obtained for the low-energy cutoff. This limit 
is always lower than the break and pivot energies and locates in the tail of the 
Maxwellian distribution of the thermal component. 
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Subject headings: Sun: flares — Acceleration of particles 



INTRODUCTION 



High-energy observations of solar flares wi th the Reuven Ramaty High-Energy Solar 
Spectroscopic Imager (RHESSI) (ILin et al.ll2002l ) allows high resolution studies over a broad 
energy range from 3 keV soft X-rays to 7-rays up to 17 MeV. The photon flux in the energy 
range of ~ 20 — 100 keV can be reasonably well fitted with a power-law function, and it s 



time-variability increases with the photon energy ( lAschwanden 



20051 : iMcAteer et all 120071 ). 



It is commonly assumed that this emission is produced by an electron population distinct 
from electrons forming a the rmal background p lasma, which is presumed to produce the low- 
energy X-ray emission (e.g. lAschwandenl |2002| ) . The impulsive high-energy emission origi- 
nates predominantly from the chromospheric footpoints, while — at least later in the fiare — 
the more slowly- varying low-energy emission is domina ted by a hot coronal so urce, observed 
in many cases to be located near EUV fiare loops (e.g. iGallagher et al.ll2002[ ). These obser- 
vations are usually interpreted in the framework of the standard fiare model where the hard 
X-ray (HXR) emission at the chromospheric footpoints of magnetic loops is bremsstrahlung 
of non-thermal high-energy electroiis mov ing downward along fiare loops from acceleration 
sites higher up in the corona ( lBrownlll97ll ). with the resulting footp oint heating and evapora- 
tion leading to the hot (usually dense) coronal thermal component ( lNeupertlll968l : iPetrosian 
I973I : lFisheiill989l ). Note, we do not automatically adopt this assumed relationship between 
accelerated and heated particles. In fact, in Section HI we interpret the observations in a 
framework where the non-thermal electrons are accelerated out of a heated thermal back- 
ground. The 'non-thermal' electron distribution is usually assumed to have a low-energy 
cutoff, the presence of which ensures that the total electron number and power are finite. 
However, it is not clear that there is a theoretical mechanism for particle acceleration which 
can naturally lead to the low-energy cutoff distinguishing non-thermal from thermal parti- 



cles (iBenz 



19771 : iMiller et al.lll997l : IPetrosian &: Liull2004l ). Indeed, it has been argued by 



Emslid (|2003[ ) that the low energy cutoff may be a redundant concept. iHannah et al.l (|2009[ ) 
suggested that a sharp cutoff in the injected electron spectrum disappears with the inclusion 
of wave-particle interactions. A dip in the electron distribution obtained through the inver- 
sion of the observed photon spectrum of some fiares may be associated with the low-energy 
cutoff. iKontar et al.l ( l2008l ). however, showed that such a feature vanishes when isotropic 
albedo correction is applied. 

The time correlation between the impulsive HXR and/or radio emission and th e deriva- 
tive of the gradual emissions at certain energies, the so-called 'Neupert effect', (INeupert 
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19681 : iDennis fc Zarrdll993l ) carries with it the imphcation that most gradual emissions are 
a "by-product", resulting from energy deposition by non-thermal electrons. This is also 
suggested by the high non-thermal electron energy content resulting from application of the 



standard collisional thick-target model (e.g. lEmslie et al.ll2004j . |2005| ). which points to high- 



energy electrons forming a dominant channel in the energy conversion process. However, 
more quantitative examination of relevant observations show that the picture is somewhat 
less clear. Some flares in volve heating of therm al coronal plasma in the absence of a power- 
law emission component (IBattaglia et al.ll2009l ). many show foot points with impulsive phas e 
emission within the energy range usually considered as thermal (iMrozek fc TomczakI 12004 ) . 
and the 'Neupert effect', which nev er represents a perfect correlation, does not hold in all 
flares or at all (thermal) energies (IMcTiernan et al.l Il999t IVeronig et al.l |2002| . |2005[ ) . So 
the possibility of heating of the solar plasma as a direct part of th e energy release before 



and /or during the ac celeration is still an open issue for investigation (jPetrosian &: Liul 12004 



Liu fc Fletcherl 120091 ). By deri ving abundance s of el ements with low first ionization poten- 
tials, such as calcium and iron, iFeldman et al.l ( l2004j ) found that at least the hot plasmas of 
some flares result from direct in situ heating of corona plasma, possibly due to a compression 
process. This approach may also lead to a measurement of the partition of hot flare plasmas 
originated from the corona and chromosphere. 

The general question of how the pre-flare magnetic energy is converted into radiation, 
plasma bulk motion , thermal and non-thermal particle energy may not have a simple answer 
(lEmslie et al.ll2005l ). Although flares share the same kind of energy source, different flares 
can have quite different appearances, and possibly involve different physical processes. Nev- 
ertheless, some well-observed characteristics can still set constraints on the overall energy dis- 
sipation process. The soft-hard-soft spectral evolutio n of some HXR pulses is one of the most 



impo rtant characteristics of high-energy emissions (IKane fc Anderson 



2004) and may point to a turbulent particle acceleration mechanism ( Grigis fc Benj 2005 ). 



1970 



Grigis &: Benz 



Early analyses (IGanl 119981 ) suggested that there is a value of photon energy at which the 
non-thermal flux does not chan ge, so that the power -law pivots about this location, a pos- 
sible further model constraint. iGrigis &: Benzl (20041) sh o wed t hat there is no single 'pivot 
energy', rather there is a small range. IBattaglia fc Benzl (|2006[ ) determined that in the rise 
phase this energy may be lower than that in the decay phase. In the context of stochastic 
particle acceleration from the thermal b ackground plasma, the pivot energy sh ould evolve 
with the background plasma properties (jPetrosian fc Liull2004j : iLiu et al.ll2010l ). However, 
distin guishing between di fferent models on the basis of observations remains a challenging 
task (IGrigis fc Benzl boosh . 



Other constraints based on the evolution of HXR light curves include the observation 
that sub-second HXR pulses peak earlier in high than in low energies, consistent with a time- 
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of-flight dispersion if the electrons producing these pulses are ac celerated at some distanc e 
from the location where the bremsstrahlung radiation is produced ( Aschwanden et al1ll996a ). 
However, this does not mean that all the energetic electrons have to be associated with these 
sub-second pulses. The reverse delay in the l onger timescale (seconds) HXR pulses could 



indicate coUisional e scape from a coronal trap (lAschwanden et al.lll996bl : lAschwandenlll998 



Krucker et al. 



Bai fc Ramatvlll979h . 



20081) . but could also be a result of a more gradual acceleration process (e.g. 



In this paper we investigate the characteristics of flare emission across a range of photon 
energies, and examine the association between temporal, spatial and spectral characteristics, 
with particular interest in the region between thermal and non-thermal parts of the spectrum. 
The paper is organised as follows. We first review theoretical considerations and present a 
simple model for the flare with an isothermal and a power-law X-ray emission component 
(Section [2]). A simple RHESSI flare on 20th February 2002, with distinct gradual low-energy 
and impulsive high-energy emissions is analysed in detail (Section |3]). An overview of the 
flare is presented in Section 13.11 The semi-calibrated photon flux is then used to derive the 
rate of change of photon fluxes at different energies during a prominent HXR pulse, and 
two temporal components are identified (Section l3.2p . This is repeated in Section 1X51 but 
using a full spectral fit. The evolution of model parameters and the corresponding photon 
fluxes are used to check self-consistency of the model, and in Section [3^ we look at the pivot 
energy derived from the rate of change of the photon fluxes. In Section HI we discuss the 
implications of these results, and conclusions are drawn in Section O 



2. Thermal and Non-thermal X-ray Emission Components of Solar Flares 

The impulsive phase of most flares is characterised by a monotonically increasing flux 
of low-energy emission and a rapidly varying flux of high-energy emission. RHESSI photon 
spectra are usually fitted with an isothermal component at low energies and a power-law 
component at high energies. For the sake of simplicity, we will ignore details of the radiative 
processes and assume that the observed photon flux consists of an isothermal component 
plus a power-law component. The photon spectrum therefore is given by 

lie, t) = ItHie, t) + iMt) (e/keV)-^(*) (1) 

where e is the photon energy, Ith{£,t) and /„,(/j(t)(e/keV)~''''^*^ correspond to the thermal 
and nonthermal component, respectively. In the presence of a pivot energy eo, Inthi^o) = 
/„t/i(eo/keV)~''' is independent of time t and the variation of the nonthermal component is 
purely due to changes in the photon spectral index 7: Inthif) (e/keV)"^^*-' = /„t/i(eo)(e/eo)~'^^*''. 
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From the spectral fitting with Ith{^) determined from the full line-plus-continuum spectrum 
derived from CHIANTI (5.2) as a function of temperature T and emission measure EM, one 
obtains /(e,to) at a given time indicated by to. The transition energy between thermal 
and non-thermal emissions, where the photon fluxes produced by the corresponding electron 
populations are equal, is determined by Ith{et,to) = /„t/i(to)(et/keV)~'^(*°\ where Inth and 7 
are fitting parameters. 

The normalized time rate of change of the photon flux at a given energy e is given by 

- J(e, t) . dt t) + htk{t) (e/keV)-^« 

Rth = hh{^,t) / hh{^,t) for e<et 

Ruth = inth{t) / lath{t) - i{t) In (e/keV) for e > 

where the dot above the relevant quantities indicates the derivative with respect to time 
and we have used the fact that the thermal and nonthermal components dominate the low- 
and high-energy photon spectra, respectively, to derive the approximate expression. The 
rate of change of photon flux of the thermal and power-law component are indicated by Rth 
and Rnthi respectively. In the RHESSI energy range we consider, the thermal continuum 
spectrum always dominates at low energies and the free- free thermal bremsstrahlung emission 
is proportional to EM{t)/ {eT{tfl'^e-x.Y>[e/kBT{t)]], where EM{t) = nl^y and V and rith 
indicate the source volume and density of the thermal electrons respectively. k-Q is the 
Boltzmann constant. 

_ ith{e,t) _ EM{t) ef f 
Ith{e,t) EM{t)^ k^T^ 2T- 

Inclusion of emission lines and free-bound emission will introduce correction terms to this 
equation. 

In principle, one can obtain the right-hand side of Equation (|2]) from the spectral fits. 
The left-hand side can be obtained from light curves of different energy bands directly. If the 
photon spectral model given by Equation ([T]) is sufficient. Equation ([2]) should be satisfied. 
If the change in the temperature of the thermal component is small, then the rate of change 
of the photon fiux is independent of the photon energy e in the thermally-dominated energy 
range (Eq. H]). One can therefore obtain Rth by fitting the rate of change of photon fiux 
in the thermal regime with a function independent of e. There will be an energy e^, where 
the time derivatives of thermal and non-thermal components are equal, i.e. Rth = Rnthi^t) 
and if the simple equivalence between non-thermal/thermal emission and impulsive/gradual 
emission holds then e[ should be comparable to obtained from spectral fits. Equation 
(2) also provides a means to investigate the presence of a constant 'pivot' energy eo for the 
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power-law component, which will occur where Rnth goes to zero - the invariant point in the 
photon spectrum of the power-law component. 



Observations 



RHESSI observed a flare on 20th February 2002 in the NOAA active region 9825, lo- 
cated near the northwest limb of the Sun at N16W80 (919"W, 285"N) a few days after its 
successful launch on 5 February 2002. It is the focus of several earlier studies d ealing with the 
characteristics during the HXR peak ( ISui et al.ll2002l : lAschwanden et al.ll2002l ). The imaging 
and spectroscopic software has been improved significantly since then and the instrumen- 
tal response is better understood and incorporated in the RHESSI software packages. We 
choose this event because of the very simple shape of its light curves, and present an inves- 
tigation here emphasizing the relationship between the low-energy and high-energy emission 
components. 



3.1. Light curves, Images, and Spectra 



This GOES 1-8 A C7.5 flare exhibits a prominent HXR pulse lasting for about 50 
seconds with count rates above 12 keV peaking near 11:06:20 UT. The first three panels of 
Figure [1] show a summary of GOES and RHESSI observations of this flare. The 3-6 keV light 
curve shows some impulsive behavior, while the 6-12 keV counts are relatively smooth. The 
RHESSI attenuator state during this flare was Al. Counts in the 3-6 keV channel when the 
attenuators are in plac e are almost all fro m higher energy photons above 11 keV, because of 
the effect of K-escape (jSmith et al.l 120021 ). We therefore only analyze counts above 6 keV in 
this work. The attenuator also reduces the count rate at low energies significantly leading 
to a livetime better than 93% for all the detectors. Pulse pileup can then be ignored in 



the spectral study (ISmith et al.l 120021 ). The rise of the RHESSI count rates below 25 keV 
becomes evident after 11:04 UT marking the onset of the flare. The slowly rising count 
rates before 11:04 UT are likely caused by particle events as is evident from the gradual 
varying count rates in higher energy channels where the statistical errors are significant. 
The background profile can be subtracted with sufficient accuracy by modehng this gradual 
varying component in different energy channels separately. 

The fourth and fifth panels of Figure [1] show the GOES temperature and emission 
measure, respectively. The background fiuxes in the two energy channels are chosen as a 
linear interpolation between average fiuxes during two intervals before and after the fiare and 
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11:00 



11:05 



11:10 11:15 11:20 

Start Time {20-Feb-02 1 1 :00:00) 



11:25 



Fig. 1. — Summary of RHESSI and GOES observations. The high-lightened range shows 
the time interval (11:05:50 — 11:06:50) used for images and spectra in Figures [2] and El 
The top panel shows the soft X-ray (SXR) lightcurves and chosen backgrounds (as linear 
interpolation of fluxes of two intervals before and after the flare) for both energy bands 
observed by GOES. The second panel shows the RHESSI count rate in several energy bands. 
The third panel shows GOES lightcurves with the background subtracted. The fourth and 
fifth panels show, respectively, the temperature and emission measure derived from the 
GOES fluxes, in agreement with results obtained by Sui et al. (2002). The sixth panel 
shows the time-derivatives of the GOES lightcurves, which is qualitatively correlated with 
the HXR count rates in the top panel. The bottom panel shows the energies radiated by the 
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are shown in the first panel. The temperature and emission measure can not be obtained 
before 11:05 UT, presumably due to the relatively low background subtracted fluxes shown 
in the third panel. The values obtained for the impulsive phase are rather i nsensitive to the 



background selection. These results are in agreement with those obtained by lSui et al.l (|2002[ ) 
for the prominent HXR pulse. We note that the temperature does not change significantly 
throughout the rise phase of the 6-12 keV count rate from 11:05 to 11:07 UT. The emission 
measure appears to grow exponentially at the beginning with a growth time of ~ 40 seconds, 
and the growth rate decreases significantly after the HXR pulse. The sixth panel shows the 
time derivatives of the GOES fluxes. Although these derivatives peak near the peak of the 
HXR pulse, and a secondary peak before the major peak appears to be correlated with 
RHESSI light curves below 25 keV in rough agreement with the Neupert effect, a broader 
correlation is not very obvious. Since the first peak does not contain high energy (> 25 keV) 
emissions which normally have a longer decay time than lower energies, we assume that 
the two peaks are independent and the first peak does not significantly affect the spectral 
properties determined for the major one. The bottom panel shows the radiative energy 
produced by the isothermal source, obt ained by fitting t he GOES fluxes using CHIANTI 



6.0.1 and assuming coronal abundances ( iDere et al.l 120091 ). The total radiated energy from 
the hot plasma is about 10^^ ergs for this flare. 

The right panel in Figure |2] shows the source structure at several energies during the 
prominent HXR peak from 11:06:10 to 11:06:30 UT obtained using the Pixon algorithm 



( iPina fc Puetterl Il993[ ) . The weak HXR coro nal source near t he solar limb has bee n inter - 



preted as the site of particle acceleration by ISui et al.l (|2002[ ). lAschwanden et al.l ( 120021 ) 



on the other hand, inferred a much smaller loop. Our results suggest that both large and 
small loops are present at energies up to 25 keV. There are also clear footpoints at > 25 keV 
which appear to be at the end of the lower energy loops. This imaging does not support 
the suggestion by Sui et al. (2002) that there is an independent low-energy thermal source 
between the two HXR footpoints, rather it indicates that the small loop connects to the 
HXR footpoints. The structure just before the flare peak is shown in the left panel of Figure 

m 

The time evolution of X-ray sources at 6-9 keV, 9-12 keV, and 12-25 keV are shown in 
Figure |3i The 6-9 keV source structure (in the left panel) is relatively simple, loop-like and 
compact before the HXR pulse. It becomes more extended in the following two time intervals 
and develops three sub-sources (though there is a possibility that this is over-resolution by the 
Pixon algorithm). The 9-12 keV emission is primarily loop- like but more extended than its 6- 
9 keV counterparts. Features associated with the footpoints start to emerge. The distinction 
between footpoint sources and loop source(s) becomes very ambiguous at 12-25 keV (in the 
right panel of Figure [3]). The brightest locations are associated with the footpoints at the 




Fig. 2. — Images of the flare for a 20 second interval before (Left) and a 20 second interval 
during (right) the HXR peak. The map color is in the unit of photons cm~^ s^^ arcsec^^. 
Left: The image and dashed contours (5, 20, and 60 % of the peak brightness) are for 9-12 
keV energy band. The solid (5, 20, and 60 %) and dash-dot (5, 20, and 60 %) contours are 
for the 6-9 keV and 12-25 keV energy band, respectively. Right: The image and dot-dash 
contours (5, 20, and 60 %) are for the 12-25 keV energy band. The solid (10, 50, and 90 %), 
dashed (10, 50, and 90%) and long-dash (2, 10, 50%) contours are for the 6-9 keV, 9-12 keV 
and 25-100 keV energy band, respectively. 
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Fig. 3. — Evolution of the X-ray images at 6-9 keV (left), 9-12 keV (middle), and 12-25 keV 
(right). The map color is in the unit of photons cm~^ s~^ arcsec"^. The image and dashed 
contours (10, 50%) are for the interval of the HXR peak. The solid contours and dash-dot 
contours are for the interval before and after the HXR peak, respectively. 

HXR peak, but the 12-25 keV structure is elongated (as is the 9-12 keV source in the middle 
panel) and may include a loop component. The X-ray images therefore reveal a complicated 
pattern of low and high energy sources, with no clear distinction between footpoints and 
loop at energies from ~10 to 25 keV. 

Figure S] shows the photon spectra fitted with an isothermal plus a power-law compo- 
nent for the three 20-second intervals of Figure [3l For this preliminary study, the RHESSI 
background is chosen as a linear interpolation, made between the average counts during the 
intervals 11:01:58 - 11:02:38 and 11:20:38 - 11:21:18. More detailed modeling of the back- 
ground is carried out in Section 13. 3^ where the spectra for 4-second intervals are analyzed. 
(Although we choose a broken power-law model in the fitting, the break energy is fixed at 
5 keV which is below the energy range of the data.) The soft-hard-soft spectral evolution 
is evident with the photon spectral index varying from 4.1 to 3.4 and to 5.5. The emission 
measure and temperature are EM = 4.7 x 10'^'^ cm~^, kT =1.5 keV; EM = 9.9 x lO^^cm"^, 
kT =1.4 keV, and EM = 1.6 x 10"^^ cm~^, kT =1.3 keV, respectively. The emission measure 
is slightly lower and the temperature is slightly higher than those from the GOES spectral 
fit, which may be attributed to the different energy ranges covered by these two instruments. 
Comparison of these spectra shows that the 12-25 keV emission is more and more dominated 
by the power-law component as the flare evolves. The spectral flt of the last interval has also 
the highest values of reduced (1-37) and residuals, leading to a probability of 8% to get 
a larger assuming a correct model. Indeed, for the last interval a thermal plus a broken 
power-law model gives much improved spectral flt. We also study the spectral evolution 
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Fig. 4. — Photon spectra of the three 20-second intervals in Figure |3] fitted with an isothermal 
plus a power-law model. Time starts from 11:05:50 and increases from left to right. Model 
parameters are indicated in the figure. 



after the impulsive phase. The simple thermal plus a single power-law model can be ruled 
out by the relatively softer high-energy spectra, which imply dominance of very low-energy 
emission by the power-law component. A thermal plus a broken power law or multi-thermal 
model can give acceptable fit. However, one should note that the high-energy photon fluxes 
change dramatically during the two intervals before and after the HXR peak and the large 
systematic residuals around the iron- line complex at 6.7 keV indicate that this feature has 
not been modeled properly. Currently, the only way to improve the modeling of emission 
lines is to fit spectra from individual detectors and take into account small gain changes and 
pulse pile up, which is beyond the scope of the current investigation. 



3.2. Semi-calibrated light curves and their rate of change 

We here study the photon flux change rate R{e, t) given by Equation ^ to quantify the 
rate of change of photon flux — how impulsive or gradual the event IS cLS cL function of 
both energy and time. We do this first in a model-independent way using the semi-calibrated 
photon flux. Semi-calibrated photon fluxes can be obtained from the observed count rates 
by using the diagonal elements of the spectral response matrix. Although the photon flux 
obtained this way does not take into account the full spectral response matrix of RHESSI, 
it can be readily obtained and gives an approximate description of the photon flux from 
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Fig. 5. — Semi-calibrated photon flux f{€,t) of tlie flare averaged over 4 s intervals. The 
curves with different colors represent photon counts in 40 different energy bands from 6 to 
50 keV. The vertical dash lines represent the start (11:05:56) and end time (11:06:44) of the 
HXR pulse. Notice that not all of the 40 energy bands are shown. 



- 13 - 



1000 ^ 



100 ^ 



1 ^ 



0.4 - 



0.2 



0.0 -r 



-0.2 - 



-0.4 



6.00 - 7.04 keV 

10.17 - 1 1.22 keV 

14.38 - 15.44 keV 

18.62 - 19.69 keV 
22.90 - 23.98 keV 
27.24 - 28.33 keV 

31.63 - 32.73 keV 
36.09 - 37.21 keV 
40.62 - 41.77 keV 
45.26 - 46.43 keV 
48.80 - 50.00 keV 



1 1 :06:00 



1 1:06:10 

Start Time 



1 1 :06:20 

(20-Feb-Oi 



1 1 :06:30 
1 1 :05:52) 



1 1 :06:40 



Fig. 6. — Top: Semi-calibrated photon flux f{e,t) and its ±lcr statistical errors during the 
HXR pulse (11:05:56 — 11:06:44). Bottom: The rate of change of the photon flux R{e,t) = 
df{e,t)/dt/f{e,t) (s~^) and its corresponding ±la errors calculated with the regularized 
method. Notice that we have shown only 11 out of the 40 energy bands. 



- 14 - 



1 .5 



0.5 



0.0 



-0.5 



-1.0 




CI = 0.05 e;= 


12.02 


C2 = 


0.29 


CI = 0.01 e;= 


10.78 


C2 = 


0.20 


C1= 0.02 £[= 


8.77 


C2 = 


0.19 


CI = 0.04 e;= 


8.80 


C2 = 


0.12 


fitting foiled :( 








C1= 0.03 e;= 


19.89 


C2 = 


-0.08 


fitting foiled :( 








CI =-0.02 e;= 


= 11.17 


C2 = 


-0.21 


C1= 0.02 e;= 


10.20 


C2 = 


-0.31 


CI = 0.01 £;= 


13.82 


C2 = 


-0.73 


C1= 0.03 c;= 


1 1.10 


C2 = 


-0.39 


CI =-0.01 e: = 


= 10.1 1 


C2 = 


-0,26 





20- 


Feb 


-02 1 1 


05 


56 


000- 


-20 


-Feb 


-02 1 1 


06 


00 


000 


20- 


Feb 


-02 1 1 


06 


00 


000- 


-20 


-Feb 


-02 1 1 


06 


04 


000 


20- 


Feb 


-02 1 1 


06 


04 


000- 


-20 


-Feb 


-02 1 1 


06 


08 


000 


20- 


Feb 


-02 1 1 


06 


08 


000- 


-20 


-Feb 


-02 1 1 


06 


12 


000 


20- 


Feb 


-02 1 1 


06 


12 


000- 


-20 


-Feb 


-02 11 


06 


16 


000 


20- 


Feb 


-02 1 1 


06 


16 


000- 


-20 


-Feb 


-02 1 1 


06 


20 


000 


20- 


Feb 


-02 1 1 


06 


20 


000- 


-20 


-Feb 


-02 1 1 


06 


24 


000 


20- 


Feb 


-02 1 1 


06 


24 


000- 


-20 


-Feb 


-02 1 1 


06 


28 


000 


20- 


Feb 


-02 1 1 


06 


28 


000- 


-20 


-Feb 


-02 1 1 


06 


32 


000 


20- 


Feb 


-02 1 1 


06 


22 


000- 


-20 


-Feb 


-02 1 1 


06 


36 


000 


20- 


Feb 


-02 1 1 


06 


36 


000- 


-20 


-Feb 


-02 1 1 


06 


40 


000 


20- 


Feb 


-02 1 1 


06 


40 


000- 


-20 


-Feb 


-02 1 1 


06 


44 


000 



10 



photon energy (keV) 



Fig. 7. — Rate of change of semi-calibrated photon flux R{e,t) = df{e,t)/dt/f{e,t) (s^^) 
versus the photon energy for 12 4-second time bins. The error bars give ±lcr uncertainties. 
The data are fitted with the model described by Equation ([5]) and the best fit models are 
shown as solid lines. For illustrative purpose, the lines are shifted vertically by the values of 
[0.6, 0.5, 0.4, 0.3, 0.2, 0.1, -0.1, -0.2, -0.3, -0.4, -0.5, -0.6] in time sequence. The fit does not 
converge for the fifth and seventh time intervals. The model parameters CI, et and C2 of 
the fit are indicated in the legend. 
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the source, and one which does not depend on assuming a particular form of the photon 
spectrum (note, we investigate full spectral fitting in Section |33|) . Following the arguments 
of Section [21 with the photon spectrogram obtained this way, one can test whether the two 
components identified from spectral fits are compatible with the two temporal components. 

Figure [5] shows some of the detailed semi-calibrated photon fluxes /(e, t) in forty energy 
bands between 6 and 50 keV. Each energy bin is set to be no smaller than 1 keV which is the 



energy resolution of RHESSI (jSmith et al.ll2002l ). The vertical axis indicates the photon flux 



averaged over a 4 second interval in the corresponding energy band. The background fluxes 
of high energy bands decrease gradually with time. We model these background fluxes for 
different energy bands with a first order polynomial fitting of background values obtained 
before and after the HXR pulse. 



We use the regularized method developed by iKontar &: MacKinnon! (|2005[ ) to obtain 
the time derivatives of these light curves. This method gives smoother derivatives while 
avoiding large errors, typical of finite differences of discrete numerical data. We assume 
statistical error for the photon flux cr/(e,t) = M^^(e) C(e, t)^^^ where M is the diagonal 
components of the instrument response matrix and C is the count rate. The top panel in 
Figure Ogives the photon flux f{e,t) and its la statistical error af{e,t) in different energy 
bands during the HXR pulse. The bottom panel shows the rate of change of photon flux 
R{e,t) = df{e,t)/dt/f{e,t) and its Icr error (TR{e,t), which is the standard deviation of 
R{e,t) modeled with a Gaussian by sampling 5000 points in df{e,t)/dt and f{e,t) within 
their respective la range of a Gaussian distribution. Here only the la range of a Gaussian 
distribution is sampled for the following reason. The relatively low flux / and its relatively 
large la error mean that a sampling over the full distribution will have points with f{e,t) 
close to zero. This will lead to very high values of R{e, t), whose distribution is poorly fitted 
with a Gaussian. The la error obtained this way should be considered as a lower limit. The 
photon fluxes at lower energies have a more gradual temporal evolution and lower rate of 
change, in contrast to the high-energy band photon fluxes, which have rapid rise and decay 
phases and highly variable rate of change. The absolute value of the rate of change often 
increases with the photon energy. However, the la error of the rate of change of photon 
flux also increases with energy, and the variation in the rates of change at high energies may 
not be significantly different from those at low energies. At the time bin for the peak from 
11:06:20 to 11:06:24, however, the rates for all energies are around zero. 

Figure [7| shows the energy dependence of the rate of change of photon flux at different 
4 second time intervals of the HXR pulse. It is clear that the higher energy fluxes have 
higher values of the derivative in the rise phase (before 11:06:20) and lower values of the 
derivative in the decay phase (after 11:06:24). This confirms that higher energy fluxes are 
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more variable than those at low energies. There appears to be two temporal components. 
At low energies, the rate of change is nearly independent of the photon energy. At high 
energies, the energy dependence of the rate of change appears to increase linearly with the 
logarithm of the photon energy. To quantify these results, we adopt the following model for 
the rate of change R 

1 Cl + C21ogio(e/e;) for e > e',. ^ ^ 

There are therefore three model parameters: CI, C2, and ej. The solid lines in Figure [7] 
indicate the best-fit model (using the curvefit function in IDL). Note that the 4th and 7th 
time intervals are not fitted with the model due to significant uncertainty in because of 
low values of the rate of change of photon flux at all energies near the HXR peak. 

A comparison of equations ([3]) and ([5]) shows that 

CI = RtH = E'M{t)/EM{t), (6) 
C2 = -lnlO-7, (7) 

where we have assumed that the temperature of the thermal component does not change 
during the HXR pulse (demonstrated using spectral fitting in the next Section), is the 
transition energy, where the rate of change of the high and low energy components are equal: 
Rth = Rnthi^t)- If the isothermal and power-law model of equation ([T]) indeed gives sufficient 
description of the observations, should be comparable to the transition energy identified 
from the spectral fit et- 



3.3. Spectral fit and the rate of change of photon flux and model parameters 

The semi-calibrated photon fiux, which is simple and fast to obtain, may give a sufficient 
approximation of the photon fiux from the source at high photon energies, but the non- 
diagonal elements of the response matrix become important at low energies. By carrying 
out a full spectral fit, which is much more time-consuming, and determining the model 
parameters from this, we can check for consistency with the results based on the semi- 
calibrated fiux. With the spectral fitting package OSPEX, the spectrum of counts from 6 
to 60 keV are fitted with an isothermal plus a single power-law model for time bins of 4 
seconds from 11:05:56 to 11:06:44 UT. Note that the model of CHIANTI (5.2) rather than 
CHIANTI (6.0.1) used in GOES data is applied here due to the lack of implementation 
of the newer model in OSPEX. Due to the change of background fiuxes in both time and 
energy, the background is separately selected for five different energy bands(3 to 6 keV, 6 to 
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Fig. 8. — Parameters of an isothermal plus a single power-law spectral fits with temperature 
adjustable (thick black lines) and temperature fixed at 1.5 keV (thin red lines). The top-left, 
top- right, bottom-left, bottom-right panels show the emission measures, temperatures, the 
power-law spectral indexes, and the of both fits, respectively. 
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Fig. 9. — Top panel: photon flux spectrum for the peak time interval: 11:06:20 to 11:06:24. 
The magenta line shows the background and the black line is the photon data with back- 
ground subtracted. The spectrum of data from 6 to 60 keV is fitted with an isothermal plus 
a single power-law model. The green line represents the thermal model and the yellow line 
is for the power-law model. The red line is for the total thermal plus power-law spectrum. 
Bottom panel: normalized residuals of the spectral fit. 
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12 keV, 12 to 25 keV, 25 to 50 keV and 50 to 100 keV) both before and after the flare and 
fitted with the first order polynomial. Since we are mostly interested in a relatively low 
energy range, where the transition between the high and low energy component occurs, an 
upper en ergy bound of 6 keV is chosen to avoid potential spectral steepening at even higher 



energies ( ISui et al. 



I2OO21) and to ensure ade quate counts above background throughout the 



Aschwanden et al.l (120021 ) also fitted the spectra of this flare with a single 



period of interest, 
power-law from 15 to 50 keV. 

We first fit with the temperature, emission measure, power-law index and normalization 
of the power-law component as free parameters. The results are indicated by the black lines in 
Figure El Both the emission measure and power-law index show significant variation during 
the HXR pulse. The variation of the temperature, however, is relatively small, between 1.3 
and 1.7 keV. To facilitate comparison with the theoretical model, we then fix the temperature 
at a typical value of 1.5 keV and do the spectral fit again. The results are indicated by the red 
lines in Figure El From the of the bottom right panel, we conclude that this model gives a 
fit to the observations which is as good as that in the model having the temperature as a free 
parameter. The smaller number of free parameters also gives smaller uncertainties for other 
model parameters, and the emission measure (top-left panel) has a smoother evolution when 
the temperature is fixed. The power-law spectral indices (bottom-left panel) from both fits 
are almost identical and both have a soft-hard-soft evolution. The goodness of the spectral 
fitting is evaluated by the reduced which is shown in the bottom-right panel. The 
for both fits are very similar except during one rise time bin (11:06:04 to 11:06:08) and one 
decay time bin (11:06:32 to 11:06:36). 

For the sake of simplicity, in the following we only use results obtained with the model 
temperature fixed at 1.5 keV. Figure HI shows the photon spectral fit with the thermal plus 
power-law model for the peak time bin (11:06:20 to 11:06:24). The emission measure is 
7.69 ± 0.18 X lO^^cm^^. The normalization of the power-law component at 50 keV is 0.66 ± 
0.01 photons s~^cm"^keV~^. The power-law spectral index is 3.3 ± 0.02. The normalized 
residuals are shown below the spectrum. The residuals are between -3 and 3 with slightly 
larger values at 6-7 keV where the iron emission lines locate. To model this feature correctly, 
one needs to fit spectra from individual detectors and take into account small gain changes 



^The time intervals chosen for fitting the background fluxes are 11:01:52 - 11:02:56 and 11:21:00 - 11:21:56 
for 3 - 6 keV, 11:02:00 - 11:02:56 and 11:19:56 - 11:21:00 for 6 - 12 keV, 11:01:56 - 11:02:32 and 11:18:04 - 
11:19:00 for 12 - 25 keV, 11:04:08 - 11:04:24 and 11:11:12 - 11:13:16 for 25 - 50 keV, and 11:04:20 - 11:04:48 
and 11:08:04 - 11:09:04 for 50 - 100 keV. We also modeled the background fluxes for eleven energy bins from 
3 to 100 keV using a third order polynomial fit with six intervals chosen around the peak time and obtained 
very similar spectral results. 
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and pulse pile up. We focus on the transition between thermal and nonthermal components 
here and leave this caveat for a future investigation. 

With the forward-fitted photon fluxes / obtained above, we calculated the normalized 
rate of change R{e,t) = {df {e,t) /dt) / f {e,t) with the same regularized method as was used 
on the semi-calibrated photon fluxes. Figure [10] shows the energy dependence of the rate of 
change with different colors representing different time bins. The errors are large for high 
energies because the photon fluxes at high energies are low and their relative errors are big. 
The solid lines are the model fit shown in Eq. Notice that the 5th and 7th lines are 
again not fitted with the model. 

According to Eqs. |6]and[7l the rate of change of the EM{t) and 7(t) can be obtained 
directly from the semi- calibrated photon flux, and fr om the spectral fits with the regularized 



method for derivatives (IKontar &: MacKinnonll2005l ). Figure [TT] compares the rate of change 



of the emission measure and the power-law spectral index with the above three methods: 
first with the semi-calibrated data (see Figure [7]), second with the photon flux derived from 
spectral fit (see Figure ITU]) , and third directly from the evolution of these parameters deter- 
mined in the spectral fit (see Figure [H]). One can see that these rates obtained with different 
methods are consistent. The rate of change of the emission measure is positive in most of 
the time bins. This indicates that the thermal emission is increasing nearly monotonically 
during the HXR pulse. The time evolution of 7 is also consistent with the soft-hard-soft 
evolution. The regularized method provides a powerful means to quantify this behavior. 



3.4. Transition and Pivot energy 

With the results obtained above, one can check the consistency of the model given by 
Equation ([1]) and test whether there is one pivot energy eo for the duration of the whole HXR 
pulse. The transition energy between the two temporal components and between the two 
spectral components should be comparable if the two spectral components have distinct 
energy and time dependences, as expected. For results obtained with the spectral fitting 
in Section |33| good spectral fits will guarantee that be comparable to et. However, with 
the semi-calibrated data, these two transition energies characterized the spectral evolution 
in two distinct dimensions. There is no guarantee that the two spectral components will 
match the two temporal components self-consistently. The top panel of Figure [12] shows et 
and obtained from the analyses above, et has a relative error of a few percent and varies 
between 10 and 12 keV, which is consistent with a constant value of ~ 11 keV. Although e't 
varies in a larger energy range of 8 to 20 keV, the relative errors are greater than 15% and 
its values are also consistent with a constant of ~ 11 keV. The agreement of these quantities 
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Fig. 10. — Same as Figured {R{e,t) in s ^) but for the photon flux derived from spectral 
modelling. 
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Fig. 11. — Top panel: comparison of the rate of change of low-energy photon flux Cli 
(red diamonds) derived with the semi-calibrated data, CI2 (blue squares) with the modelled 
spectral data, and the rate of change of the emission measure EM /EM (green circles). 
Bottom panel: comparison of parameters C2i (black triangles) and C22 (orange crosses) 
derived with the semi-calibrated data and modelled spectral data respectively. The rate of 
change of the spectral index 7 can be obtained from C2i and C22 as shown in the legend. 
Red diamonds, blue squares and green circles represent 7 from semi-calibrated photon fluxes, 
from modelled spectral photon fluxes, and from the time derivatives of the power-law index 
respectively. Error bars of EM /EM and d'j/dt (both in green circles) are obtained with the 
regularized method. Error bars of all the other parameters indicate the la uncertainties of 
curve fit in Figures [7] and [TH] and described by Eq. . 
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Fig. 12. — Top panel: evolution of the transition energy between the low and high energy 
component derived from variability of the semi-calibrated data (red diamonds) , the modelled 
spectral data (blue squares) and from the isothermal plus power-law spectral fit (green 
circles). The error bars of e^^ and e^2 show, respectively, the uncertainties of curve fit in 
Figures 171 and fTOl as described by Eq. ([5]). The error of et is the standard deviation of 5000 
simulated intersections between thermal and nonthermal components taking into account 
uncertainties of all the spectral fitting parameters. Also shown is the possible minimum 
electron low-energy cutoff Ec (Section H]). Bottom panel: pivot energy cq with la errors (see 
Section 13.41 for more details) derived from the variability of the semi-calibrated data (red 
diamonds) and modelled spectral photon fluxes (blue squares). 
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implies consistency of the model. The big error bars of e[ are due to the uncertainty of 
determining the cross point of the two lines in the fitting shown in Figures [7| and [TOl This 
may reveal a complicated physical process where the transition between the slow-varying 
gradual component and the impulsive component is rather an energy range than a single 
point. 

The usual method of determining eo by spectroscopic fitting depends on the assumed 
spectral model. We instead use the above rate of change study to derive eo and its errors for 
each interval. From equation ([S]), one can show that the rate of change of the photon flux of 
the power-law component is zero at 



This is the pivot energy at a given time interval. With parameters CI, C2, and obtained 
above, we calculated the pivot energy eo as shown in the bottom panel of Figure [121 The 
error of eo is taken as the standard deviation of simulated eo with 5000 sampling points 
of ej, CI, and C2 distributed within their respective la range of a Gaussian distribution. 
The results are consistent with a constant value of ~ 9 keV except for the fourth and tenth 
intervals where deviations of eo from 9 ke V are greater than Za. T hese values of the pivot 
energy are also less than those obtained by iBattaglia fc Benzl (120061 ) for the loop top source 
of a few other flares, but are in agreement with these of the footpoint sources especially for 
values in the HXR decay phase. Since CI is mostly positive and C2 evolves from positive 
to negative values from the HXR rise to decay phase, lO"*^^/*^^ evolves from less than 1 
to greater than 1, implying that eo increases from < e[ to > e[. It should be noticed 
that when C2 approaches zero, the amplitude of C1/C2 approaches infinity and lO"*^^/*^^ 
approaches either or infinity, both of which are not physical. It also leads to huge error bars 
for time bins near and after the HXR peak. Indeed, a sampling of C2 over a full Gaussian 
distribution after the peak (when C2 < and CI > 0) will lead to infinite values of eo, whose 
distribution is poorly fitted with a Gaussian. Theoretically, for HXR pulses with soft-hard- 
soft spectral evolution, if the transition energy and temperature do not change significantly 
and the emission measure has a gradual and monotonic increase, the pivot energy in the 
decay phase should be higher than that in t he rise phase, as indi c ated in Figure [T2l This 
result is in agreement with previous studies flGrigis fc Benzl 120051 : IBattaglia &: Benzl |2006[ ) 
and may be attributed to an effect of chromospheric evaporation (iLiu et al.l l2010l ). The 
pivot energy eo is also comparable to the transition energies et and ej, in agreement with 
the scen ario where energetic electrons are accelerated froni a low -energy thermal background 
plasma jBenzlll977l : IPetrosian k Liul 120041 : iGrigis k Benzl booeh . 
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Discussion 



We started with the hypothesis that there are two distinct emission components with 
low energy photons evolving gradually and high energy photons having a rapid evolution. 
One consequence of this is that, as long as the temperature of the thermal component varies 
slowly (much slower than the emission measure), which simplifies the model significantly [see 
eq. dl])], the break energy between non-thermal and thermal emission in the photon spectrum 
should be comparable to the transition energy between slowly- and rapidly- varying photon 
fluxes found by evaluating time derivatives. Within the uncertainties of this method, we have 
demonstrated that this is the case (top panel of Figure [12]) , and that the transition energies 
are always around 11 keV. However, it is clear that there are substantial error bars on the 
values of the transition energy e^, which are relatively independent of whether the spectral 
fitting approach (model-dependent) or the semi-calibrated approach (model independent) is 
used. In fact it is not possible to pin down the gradual/rapid boundary within about ±5 keV 
throughout most of the flare, especially near the HXR peak. This is due to the difficulty of 
determining the folding point of the broken line. Better data with much higher count rates 
and lower statistical errors are required — for example, a more intense but equally simple 
flare — to examine whether or not such a boundary can be more clearly identified. 

To understand further the relationship between the high and low energy emission compo- 
nents produced presumably by two distinct electron populations through the bremsstrahlung 
process, we investigate the electron numbers in each population. In one version of the stan- 
dard model, it is postulated that electrons are accelerated at a reconnection current sheet and 
the acceleration process is decoupled from t he electron transport and magnetic field evolution 
after the reconnection ( Aschwanden Il998), which woul d imply an ideal Neupert effect not 
wholly supported by observa tions (IVeronig et al.l 120051) . This scenario also encounters the 
well-known number problem (Fletcher &: Hudson 20081 ). Given the high energy release dur- 



ing some large flares (lEmslie et al 



2004 



20051 ). energy flows likely play more important roles 



than nonthermal electron fluxes in our exploration of the physics in the impulsive phase. It is 
possible that a significant fraction of the magnetic energy is converted into particle energies 
after the reconnection during the relaxation of magnetic field lines. The reconnection only 
permits the changes of magnetic field topology and may not correspond to the dominant en- 
ergy dissipation and parti cle acceleration process, which can proceed after the reconnection 
([Fletcher fc HudsonI 120081 ). In the context of stochastic particle acceleration, it is usual to 
assur ae that electrons arrivin g at the chromosphere are accelerated out of a population in the 
loop (jPetrosian fc Liu! [2004 ). therefore — assuming for simplicity no magnetic convergence 
— the n umber density of non-thermal electrons should be no larger than the loop number 
density. [Sui et al.l (120021 ) suggested that high-energy electrons might be accelerated from 
a cool background plasma not observed in X-rays. There is no observational evidence for 
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such a cold background. Theoretically, it is also difficult to understand why the acceleration 
should proceed in relatively cool regions given the microphysics of particle energization by 
electric fields is the same for both thermal and nonthermal populations. Moreover, as we will 
show below, the observed thermal plasma is dense enough to provide electrons responsible 
for the high energy emission. The assumption of a cooler background source for the high 
energy electrons appears to be unnecessary. In the following, we will assume that electrons 
producing the high energy component are accelerated from the observed thermal component. 
The number density of a nonthermal electron beam can be estimated as 



^nth 



EVeAnXR ' 



(9) 



where Pc is the power in electrons of energy greater than Ec, Ahxr is the area of HXR 
footpoints where electrons enter the chromosphere and can be estimated from fiare images, 
E is the average electron energy and Ve is the corresponding electron velocity. Pc is given by 



Pr 



At 



5-2 



(E,/keV)~(^-2)[keVs 



(10) 



where 5 is the spectral index of the underlying non-thermal electron fiux spectrum (5 = 7 + 1 
in the collisional thick target r nodel), is a normalization parameter and is numerically re - 
lated to the photon spectrum (jBrownlll971uSaint-Hilaire fc Benzll2005l : iFletcher et al.ll2007l ): 



Ar 



6.44 X 10 



33 7(7 - 1) . 

5(7-1,1/2) ^' 



(11) 



where B is the beta function, and A^ is the normalization of the power-law fit to the photon 
spectrum /(e) = ^^(e/keV)"''' (in photons s~^cm~^keV"^). It can be shown that E = 
F{E)EdE/ /^^ F{E)dE = {S -l)Ej(S-2) where F{E) = AE{E/keV)-^ is the electron 
fiux injected into 'the footpoints JSrownl EoTlh . Ve = J^^F{E)dE/ j^^F{E) /ve{E)dE = 
{6 — l/2)ve{Ec)/{S — 1), where we have assumed the nonthermal electron distribution is 
given by F{E)/{ve{E)AHXR). 

The thermal electron number density is estimated as nth = {EM/VsxrY^'^ , where EM is 
the emission measure of the thermal component and Vsxr is the volume of the SXR thermal 
coronal loop. For the sake of simplicity, Vsxr is evaluated as where Asxr is the 

projected area of the observed thermal coronal source. The areas within the 30% contours 
of the maximum value of 25-100 keV (for HXR source) and 6-9 keV (for SXR source) images 
can be obtained from the Pixon images directly. From the peak-time 20-second integrated 
image shown in Figure [2l we estimate Ahxr to be 20 square arcsec, Asxr to be 100 square 
arcsec, and Vsxr is then about 1000 cubic arcsec. The corresponding rith is greater than 
~ 10^^ cm-3. 
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The fraction of the thermal electrons accelerated into a non-thermal distribution a 
should be less than 1. We then have Unth < and 



Setting a = 1, which would correspond to the minimum possible value of i?c, and with the 
parameters EM, 7, A^ from Section 13.31 we calculate the low limit of E^ in each 4 second 
time bin of the spectral fitting. A comparison of low limit of E^. and the transition energies 
derived from variability of the semi-calibrated data (e^^^), the modelled spectral data (e^g) 
and from the isothermal plus power-law spectral fit (e^) is presented in Figure [121 At the 
peak of the flare, the minimum possible E^. is around 5.2 keV, and throughout the event it is 
always smaller than the transition energy et, which is around 11 keV. This is fully consistent 
with the stochastic acceleration model where non-thermal electrons are accelerated from a 



thermal background (jBenzlll977l : iPetrosian fc Liull2004j : iGrigis fc Benzll2006l ). We note that 5 



is always greater than 4 for this flare, the E^. obtained with equation ( !T2|l is rather insensitive 
to the poorly-determined source size Ahxr and Asxr- 

Rearranging Equation f|T2|) . and keeping other factors constant, a and E^ vary as a oc 
Ec If we set the low-energy cutoff E^ ^ et ^ 10 keV, with the low limit of E^. = 5 

keV at a = 1 for the flare peak where 5 = 7 + 1 = 4.3, this gives a = 0.07. In other words, 
accelerating 7% of the hot thermal distribution in the loop would satisfy the footpoint 
requirements. One may further assume that these accelerated electrons come from the high 
energy tail of the thermal distribution. In the tail of a Maxwellian, the fraction of electrons 
with energy above E = xkT is a^^ = A/4x/vrexp(— x), so a = 0.07 corresponds to x = 3.4. 
Since kT = 1.5 keV, the non-thermal population would correspond to the accelerated tail 
of electrons with initial energy above 5.1 keV. This in turn sets a requirement that the 
acceleration timescale would have to be less than the electron-electron Coulomb collision 
timescale for an electron of energy 5.1 keV in a Maxwellian plasma of temperature 
1.5 keV and density 4.5 x lO^^cm^'^ (using EM, Vsxr determined above). The value of 
Tee for an electron in the core of such a distribution is 0.02 s, and that for an electron of 
energy x^T is approximately x^^'^^ee or 0.13 s in this case. Note that this threshold energy 
should not be compared with the transition energies in the overall photon spectra directly. 
The electron transport and X-ray emission processes will make the transition energy in the 
electron population different from that of the emitted photons. In the context of stochastic 
acceleration in the flare loop, a higher break energy in the overall photon spectrum than in the 
electron distribution in the corona acceleration site implies that the electron escape timescale 
from the acceleration site to the footpoints is relatively long compared to the Coulomb 
collision ti mescale at these break energies so that the HXR fluxes from the footpoints are 



suppressed (IPetrosian &: Liull2004l ). 
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It should here be remarked that the low values obtaine d for Ec call into question the 



application of the cold coUisional thick target. lEmslid ( 120031 ) shows that the cold collisional 



thick target loss rates are overestimated for electrons of energy less than 5 kT. Though we 
have put bounds on E^, it remains a fit parameter, and one can instead consider the energetics 
of an injected non-thermal electron distribution extending from kT, which merges into the 
ambient thermal distribution. For this event at its peak (7 = 3.3) the collisional thick target 
power requirement above 5kT is 7.3 x lO^^erg s~^, which is already too high compared with 
the total radiatio n energy of ~ 10^^ ergs obtained from the GOES observation. Following 



the calculation of lEmslid (|2003[ ) gives about 35 times as much as this in total injected power 
above kT, which further demonstrates the necessity of going beyond the classical cold thick 
target model. 

The discussion above assumes a beam of electrons with a power-law distribution. Con- 
sidering the pitch-angle scattering, the average electron velocity will be lower giving rise to a 
higher local non-thermal electron density and hence a larger Ec- Indeed, X-ray images reveal 
more complex structure near the transition energy between thermal and non-thermal compo- 
nents. The separation of the emission into two distinct components is rather ambiguous near 
the transition energy. It is possible that this separation is an artifact of a simplified model of 
X-ray spectrum and the energy dependence of the rate of change of the photon fluxes. Two 
distinct electron populations can be due to different physical processes, which dominate at 
different energies, with the electron behavior varying considerably from low to high energies 
and the apparent distinction between the low and high energy emission components being 
just a consequence of these processes. 



5. Conclusions 

We have developed a new method to study in detail the temporal evolution of thermal 
and non-thermal photon fluxes in solar flares. The application of this method to a flare on 
Feburary 20 2002 demonstrates that as expected, the low energy part of the spectrum evolves 
slowly, and the high energy part evolves rapidly, with an intermediate range between a few 
keV and 20 keV where the behavior is in transition. The data support the scenario in which 
the non-thermal component of the flare spectrum is impulsive, and the thermal component 
is gradual, in that the transition energies between these two behaviors are the same within 
errors whether examined in time or in energy. 

However, although in the spectral fitting exercise it is possible to make a clean separation 
between a non-thermal, impulsive component, and a thermal, gradual component, time 
evolution gives a more ambiguous picture, due to the large error bars. Imaging is also 
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ambiguous, with no clear distinction between footpoints and loops in the energy range around 
9-25 keV. Therefore we must leave open the possibility that the electrons form a continuous 
distribution over this range. Further studies with larger flares should help to improve the 
precision with which we can identify the transition between gradual and impulsive behavior. 

Finally, the presence of a single pivot point throughout the flare is not supported by 
our analysis, though a pivot 'range' is. There is some evidence of a slightly higher value for 
this pivot range in the decay phase than in the rise phase. 
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